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1.  Introduction 


Fracture  in  engineering  materials,  particularly  metals,  often  initiates  at  inclusions  and  voids  (1). 
Figure  1  is  a  scanning  electron  micrograph  of  the  fracture  surface  from  a  failed  type  431  stainless 
steel  bolt  from  a  helicopter.  The  dark,  geometric  shapes  on  the  fracture  surface  are  grains,  and 
the  appearance  of  the  fracture  surface  is  typical  of  brittle,  intergranular  fracture.  The 
lighter-colored  amorphous  objects  on  the  fracture  surface  are  inclusions.  As  figure  1  illustrates, 
cracks  often  interact  with  inclusions  and  voids.  This  interaction  is  a  real-world  problem  and 
much  research  has  been  conducted  since  the  early  1960s  in  the  fracture  mechanics  area  of  the 
crack-inclusion  interaction.  Developing  an  understanding  of  the  crack-inclusion  interaction  is 
important  to  mechanical  design,  life  cycle  predictions,  failure  analysis,  and  many  other 
engineering  considerations.  This  report  reviews  some  of  the  foundational  work  in  the 
crack-inclusion  interaction  of  fracture  mechanics. 

There  are  many  types  of  inclusions  in  engineering  materials.  In  metals,  we  find  casting  defects, 
dislocations,  different  phases,  impurities,  and  voids.  For  example,  silicon  carbide  particles  can  be 
inclusions  in  aluminum.  Ceramics  often  contain  various  phases  such  as  graphite  inclusions  in 
diamond.  In  composites,  we  may  consider  the  reinforcing  materials  as  inclusions  within  a 
matrix,  for  example,  carbon  fibers  in  an  epoxy  matrix.  In  building  materials,  we  see  sand  and 
gravel  as  inclusions  in  cement  and  aggregate  of  various  sizes  in  concrete  and  bricks,  and  adobe  is 
full  of  sand  inclusions  within  a  clay  matrix.  Building  materials  are  also  riddled  with  voids  and 
porosity.  Polymers  contain  varying  degrees  of  crystallinity  and  varying  lengths  of  polymer 
chains,  and  these  too  may  be  considered  inclusions.  Even  biological  materials  contain 
inclusions — for  example,  voids  in  the  foam-like  structures  of  bones  and  bird  feathers  may  be 
thought  of  as  inclusions.  Perhaps  we  could  say  any  inhomogeneity  in  engineering  materials  is  an 
inclusion.  Since  we  know  there  is  no  truly  homogeneous  material,  all  materials  may  be  said  to 
have  inclusions,  so  the  importance  of  understanding  the  interaction  between  a  crack  and  an 
inclusion  becomes  clear.  But  to  simplify  our  considerations,  we  may  think  of  inclusions  as  either 
different  phases  within  a  material  or  void,  which  is  the  absence  of  material.  In  accordance  with 
this  idea,  we  consider  inclusions  to  be  divided  into  these  two  categories:  voids  and  phases. 

Figure  2  is  a  micrograph  of  a  sample  of  poorly  welded  titanium  (Ti-6A1-4V).  To  avoid  oxidation, 
the  alloy  must  be  welded  in  an  inert  atmosphere  and  so  typically  argon  is  flooded  over  the 
titanium  as  it  is  welded.  During  the  welding  process,  the  grains  of  metal  melt,  recrystallize,  and 
grow.  If  the  process  is  poorly  controlled,  voids  may  develop  as  the  inert  gas  is  trapped  within  the 
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Figure  1.  Micrograph  illustrating  interaction  between  fracture  surface  and  inclusions. 


structure  or  bubbles  up  to  the  surface.  Figure  3  shows  a  region  of  heat-treated  steel  that  has 
undergone  a  martensitic  transformation.  The  pre-heat-treated  micro  structure  of  steels  is  often 
largely  composed  of  cementite  (Fe3C)  and  ferrite  (a  phase  iron).  Heat  treatment  of  steels  causes 
phase  transformation  to  austenite  (7  phase  iron),  and  if  the  steel  is  quenched,  the  austenite  will 
transform  into  martensite.  Martensite  grains  are  much  smaller  than  austenite  grains,  and  this 
contributes  to  the  increased  strength  of  the  martensite  structure  (see  the  Hall-Petch  equation, 
which  relates  the  strength  of  a  metal  to  the  reciprocal  of  the  square  root  of  its  grain  size,  for 
example,  Reed-Hill  and  Abbaschian,  (1 ),  which  has  an  excellent  discussion  on  steel 
microstructures  and  strength).  The  martensitic  transformation  also  causes  an  increase  in  volume, 
and  this  is  why  manufactured  steel  parts  are  often  machined  to  near  net  shape  prior  to  heat 
treating  and  then  machined  to  the  final  shape  after  heat  treatment.  However,  sometimes  the 
martensitic  transformation  may  remain  incomplete  such  that  there  is  a  small  amount  of  retained 
austenite.  The  bright  white  spots  in  figure  3  are  retained  austenite.  Obviously,  these  two  phases 
of  steel  will  have  different  elastic  properties.  How  will  a  crack  behave  when  it  approaches  an 
inclusion  that  is  less  stiff  than  the  surrounding  material?  How  will  a  crack  behave  if  it 
approaches  an  inclusion  that  is  stiffer  than  the  surrounding  material? 
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Figure  2.  Micrograph  of  voids  in  titanium. 
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Figure  3.  Micrograph  of  retained  austenite  (white)  in  martensitic  steel. 
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Figure  4  shows  small  circular  carbide  inclusions  in  steel,  as  seen  in  a  cross  section  perpendicular 
to  the  primary  fracture  surface.  The  secondary  crack  on  the  left  is  approaching  a  carbide 
inclusion  and  the  secondary  crack  on  the  right  has  traveled  around  a  carbide  inclusion.  A  crack 
may  pass  near  an  inclusion,  perhaps  attracted  to  the  inclusion;  or  a  crack  may  move  away  from  an 
inclusion,  perhaps  repelled  by  the  inclusion;  or  a  crack  may  cross  the  boundary  into  and  perhaps 
propagate  through  an  inclusion.  Additionally,  we  might  consider  whether  the  inclusion  is 
perfectly  bonded  to  the  matrix,  partially  bonded  to  the  matrix,  or  not  bonded  to  the  matrix.  Each 
of  these  considerations  affects  the  result.  In  this  study,  we  assume  that  a  crack  passes  near  an 
inclusion,  and  the  inclusion  is  perfectly  bonded  to  the  surrounding  matrix. 


Figure  4.  Micrograph  of  carbide  inclusions  in  fractured  steel. 


2.  Inclusions  and  Cracks 


According  to  King  (2),  renowned  scientist  John  D.  Eshelby  was  a  “pure  theoretician”  who  would 
rather  reason  his  way  to  the  answer  than  conduct  experiments.  By  1961,  Eshelby  had  authored 
seminal  works  on  the  mechanics  of  inclusions  (5,  4).  Famous  for  his  thought  experiments, 
Eshelby  determined  the  stress  field  around  an  ellipsoidal  inclusion  using  the  superposition  of  two 
linear  elastic  problems  solved  using  a  Green’s  function  approach  (5).  He  imagined  a 
homogeneous  linear  elastic  solid  volume,  for  example,  figure  5a,  with  some  known  elastic 
constants,  say  El]ki.  Within  this  semi-infinite  volume  he  imagined  a  smaller,  ellipsoidal 
volume — the  inclusion.  The  ellipsoidal  volume  experiences  a  uniform,  inelastic  deformation,  for 
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example,  a  martensitic  transformation  as  discussed  earlier,  which  causes  the  inclusion  to  grow  in 
volume.  Since  the  inclusion  is  confined  by  the  surrounding  matrix,  the  deformation  puts  the 
inclusion  and  the  matrix  into  some  state  of  stress.  Eshelby  imagined  first  removing  the  inclusion 
from  the  matrix,  as  seen  in  figure  5b,  where  the  dotted  line  is  the  shape  of  the  inclusion  when  it  is 
confined  by  the  matrix.  Once  removed  from  that  confinement,  the  now  free  inclusion  will  relieve 
the  deformation  and  seek  equilibrium.  Then  with  no  force  on  the  matrix  volume,  Ym,  or  the 
inclusion  volume,  V/,  the  inclusion  will  assume  a  uniform  strain,  which  is  known  as  the 
eigenstrain,  e*-,  that  is,  the  strain  at  zero  stress.  Then  the  strain,  stress,  and  displacement  in  the 
matrix  and  inclusion  are  given  by  equations  1  and  2,  respectively. 


(c)  (d) 


Figure  5.  Eshelby’s  inclusion  analysis. 


Eij  0  , 

Oij  =  0  ,  (1) 

Ui  =  0  . 
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(2) 


£ij  ~  £ij  i 
@ij  0  , 

U%  &ijXj  . 

Next,  Eshelby  imagined  applying  a  traction,  Tj,  to  the  surface  of  the  inclusion  sufficient  to  force 
the  inclusion  into  its  original  shape,  the  dotted  line  in  figure  5c.  This  traction  induces  an  elastic 
strain  in  the  inclusion  that  exactly  cancels  the  eigenstrain,  that  is,  £®J  =  —  c*;.  The  traction  on  the 
inclusion  surface  is  given  by  equation  3,  where  a*;]  is  the  stress  computed  from  the  eigenstrain. 
Then  the  strain,  stress,  and  displacement  in  the  matrix  are  still  given  by  equation  1,  but  the  strain, 
stress,  and  displacement  in  the  inclusion  are  now  given  by  equation  4. 


Tj  OijTlj  aijnj 


(3) 


£  ■ .  —  -i_  c* .  —  n 

cij  '  cij  u  > 
<Jij  =  EijklSij  =  —EijklSij  =  ~  O’y  1 

Ui  =  0  . 


(4) 


With  the  inclusion  deformed  back  to  its  original  shape  by  the  application  of  traction  Tj,  Eshelby 
imagined  that  the  inclusion  was  placed  back  into  the  matrix.  There  is  no  change  in  the  strain, 
stress,  and  displacement  of  the  inclusion  since  the  traction  Tj  is  still  applied;  that  is,  the  state  of 
stress  of  the  matrix  is  still  given  by  equation  1  and  the  state  of  stress  in  the  inclusion  is  still  given 
by  equation  4.  But,  if  traction  Tj  were  removed  from  the  inclusion,  the  system  would  return  to  its 
original  strain,  stress,  and  displacement  state,  which  is  not  known.  So  instead  of  removing  the 
known  traction  Tj  that  was  applied,  Eshelby  applied  a  canceling  force  Fj  =  —Tj  to  the  internal 
surface  of  the  matrix,  the  surface  that  was  exposed  when  the  inclusion  was  removed,  surface  S /  in 
figure  5d.  Then  Eshelby  was  able  to  quantify  the  displacement  field  of  the  elastic  body,  called  the 
constrained  displacement  field,  in  terms  of  the  applied  force  Fj  on  the  surface  S /,  using  the 
Green’s  function  of  the  elastic  body,  that  is,  the  matrix. 

From  equation  3  and  Fj  =  —Tj  we  have 
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Fj  =  -Tj  =  a*knk  . 


(5) 


With  the  Green’s  function,  GtJ,  for  the  elastic  body,  the  constrained  displacement,  u-',  is  given  by 


uci(xi)  =  [  Fj(x,i)Gij(xi,x'i)dS(x'i)  ^  I  a*jknk(x'i)Gij(xi,x'i)dS(x'i)  .  (6) 

iS;  JSt 

We  know  that  strain  is  the  derivative  of  displacement  (the  displacement  gradient),  so  the  strain  in 
the  elastic  body  is  given  by 


0tknk(^Guj(xi,x'JdS(x'J  , 


(7) 


<Tiknk(x3[Guj(xi,x!J  +  Gjiii(xi,x!J]dS(x'i)  , 


(8) 


and  stress  in  the  elastic  body  due  to  the  constrained  displacement  field  is  given  by 


Vijfa)  =  Eijkl£ckl(Xi)  .  (9) 

Finally,  we  can  solve  for  the  state  of  stress  in  the  inclusion  from  the  state  of  stress  in  the 
surrounding  matrix.  The  displacement  and  strain  will  be  the  same  since  the  inclusion  is  assumed 
to  be  perfectly  bonded  to  the  surrounding  matrix.  The  stress  in  the  inclusion  will  be  equal  to  the 
stress  in  the  surrounding  matrix  minus  the  stress  in  the  inclusion  due  to  the  applied  traction. 
Therefore,  the  strain,  stress,  and  displacement  in  the  matrix  is  given  by  equation  10  and  the  strain, 
stress,  and  displacement  in  the  inclusion  is  given  by  equation  11. 


£ij  —  eij  t 
°ij  =  aij  ) 
Ui  =  Ui  . 


(10) 
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(11) 


f  ■  ■  —  . 

^  ^  =  Eijki(4i  -ell), 

Ui  =  Ui  . 


We  now  have  a  linear  elastic  solution  for  the  strain,  stress,  and  displacement  both  inside  and 
outside  of  an  elliptical  inclusion  in  an  infinite  elastic  body.  Further  mathematical  work  involving 
Eshelby’s  tensor,  Stjki,  which  relates  the  constrained  strain  to  the  eigenstrain  as  e\3  =  Sijki£*ki , 
may  be  needed  to  identify  explicit  expressions,  but  this  is  discussed  in  more  detail  by  Weinberger 
et  al.  (5).  What  we  have  learned  is  that  all  we  have  to  do  to  obtain  numerical  values  for  these 
states  of  stress  is  identify  the  Green’s  function  for  the  elastic  body. 

But  we  are  interested  in  how  a  crack  interacts  with  an  inclusion.  So  let  us  consider  how  the  stress 
field  is  affected  by  the  presence  of  a  crack.  If  we  know  the  solution  for  the  stress  field  in  an 
infinite  elastic  body  due  to  the  presence  of  a  crack,  then,  as  in  figure  6,  we  can  superpose  that 
solution  on  the  solution  we  just  found  for  the  stress  field  of  an  infinite  elastic  body  with  an 
inclusion.  Therefore,  if  we  want  to  know  how  a  crack  will  interact  with  an  inclusion,  all  we  need 
is  the  stress  solution  for  a  cracked  infinite  elastic  body  and  the  Green’s  function  for  an  infinite 
elastic  body  with  an  inclusion. 


I  1  I 


Figure  6.  Superposition  of  solutions  from  a  cracked  elastic  body  and  an  elastic  body  with  an  inclusion. 
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3.  Dislocations,  Green’s  Functions,  and  Superposition 


3.1  Dislocations  and  Inclusions 

A  dislocation  is  a  line-shaped  defect  in  a  crystal  lattice  where  there  is  mismatch  in  the  crystal 
structure.  Edge  dislocations  are  formed  by  deformation  of  a  body;  they  he  perpendicular  to  their 
Burgers  vector  and  move  in  the  direction  of  the  Burgers  vector  through  the  body  under  load.  For 
a  much  more  detailed  discussion  on  dislocations  and  Burgers  vectors,  see,  for  example,  Reed-Hill 
and  Abbaschian  (7)  or  Tadmor  and  Miller  ( 6 ). 

In  1964,  John  Dundurs  and  Toshio  Mura,  both  professors  at  Northwestern  University,  described  a 
solution  for  an  arbitrarily  oriented  single  edge  dislocation  interacting  with  an  inclusion  in  an 
elastic  body  (7).  They  found  the  solution  by  assuming  the  Airy  stress  potential  and  solving 
backwards  to  determine  the  terms  by  intuition  and  trial  and  error  (8). 

Dundurs  and  Mura  provided  the  first  solution  in  linear  elasticity  for  an  edge  dislocation  in  an 
infinite  elastic  medium  containing  a  circular  inclusion,  which  has  elastic  properties  different  from 
the  surrounding  matrix  (7).  They  assumed  a  perfect  bond  between  matrix  and  inclusion  so  that 
the  displacement  and  stress  fields  are  the  same  for  both  matrix  and  inclusion  at  the  boundary 
between  the  matrix  and  the  inclusion.  Further,  they  assumed  there  is  a  singularity  at  the  edge 
dislocation,  but  no  singularity  anywhere  else  in  the  matrix  or  inclusion.  Also,  any  closed  contour 
integral  that  does  not  contain  the  center  of  the  dislocation  is  required  to  equal  zero,  but  any  closed 
contour  integral  that  does  contain  the  center  of  the  dislocation  is  required  to  equal  the  Burgers 
vector  of  the  dislocation,  b.  Finally,  the  stress  at  infinity  is  assumed  to  vanish. 

The  complex  potential  method  of  solution  for  a  dislocation  interacting  with  an  inclusion  is 
beyond  the  scope  of  this  report.  Dundurs  and  Santare  (8)  provide  an  excellent  discussion  of  the 
topic  and  a  review  of  the  literature. 

3.2  Dislocations  and  Green’s  Functions 

Dundurs  and  Santare  (8)  provide  a  concise  discussion  on  using  a  “concentrated  unit  of 
inhomogeneity”  as  a  forcing  function  in  the  solution  to  a  boundary  value  problem,  which  they 
attribute  to  Stakgold  (9)  who  first  wrote  of  this  possibility  in  1967.  Known  as  a  Green’s  function, 
Dundurs  and  Santare  (8)  discuss  that  the  solution  to  this  concentrated  unit  of  inhomogeneity,  the 
Green’s  function,  may  be  used  to  find  the  solution  for  a  general  boundary  value  problem  in  terms 
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of  a  continuous  forcing  function.  The  boundary  value  problem  they  discuss  is  the  displacement, 
strain,  and  stress  fields  and  the  corresponding  boundary  conditions  found  in  elasticity  problems. 
Since  the  forcing  function  is  provided  by  a  concentrated  unit  of  inhomogeneity,  we  may  consider, 
for  example,  a  point  load  in  an  elasticity  problem,  or,  more  appropriate  to  this  discussion,  a 
dislocation.  Dundurs  and  Santare  provide  a  concise  mathematical  discussion  of  the  use  of  a 
dislocation  as  the  required  singularity  in  an  otherwise  undisturbed  elastic  body,  which  provides 
Green’s  function,  which  leads  to  the  stress  field  solution  of  the  boundary  value  problem  in  plane 
elasticity  (S). 

What  is  a  Green’s  function1?  To  answer  that  question,  first  we  must  understand  the  Dirac  delta 
function.  The  Dirac  delta  function  is  used  to  describe  a  sudden,  rapid  rise  in  force,  described  by  a 
forcing  function.  So,  for  example,  if  we  have  an  undisturbed  piece  of  plywood  lying  on  the  floor 
and  we  place  a  nail  somewhere  in  the  middle  but  do  not  yet  strike  the  nail  with  a  hammer,  then  the 
plywood  is  at  a  state  of  zero  stress  everywhere.  Suddenly,  we  strike  the  nail  with  a  hammer  so 
that  the  plywood  is  at  a  state  of  zero  stress  everywhere  except  under  the  tip  of  the  nail  where  the 
stress  rapidly  approaches  some  huge  value  (call  it  a  conceptual  infinity).  The  forcing  function 
describes  the  force  applied  by  the  nail  to  the  plane  of  plywood,  and  the  Dirac  delta  function 
describes  the  change  in  condition  from  zero  everywhere  to  infinity  under  the  tip  of  the  nail. 

The  Dirac  delta  function  has  these  properties: 


5(x-0 

5  (x  —  £)  dx 


0, 

1, 


f(x)S  (x 


0  dx  =  /(£)> 


x  ±  i , 

/3>0, 

P  >  0. 


(12) 

(13) 

(14) 


Mathematically,  to  solve  a  differential  equation  containing  the  Dirac  delta  function,  we  need  to 
transform  the  Dirac  delta  function.  We  use  the  third  property,  equation  14,  to  perform  a 
transformation  using  a  linear  differential  operator,  C,  so  that 

'Apparently,  there  is  some  debate  as  to  whether  it  is  correct  to  say  “Green  function”  or  “Green’s  function”.  A 
somewhat  recent  article  showed  that  the  occurrence  of  the  use  of  “Green’s  function"  has  taken  precedence  in  recent 
years  (JO).  Since  Green’s  functions  are  problem  specific  and  more  idea  than  identity,  I  prefer  to  use  “Green’s  func¬ 
tion". 
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£9{x,€)  =  S{x-£)  ■ 


(15) 


According  to  Stakgold  (9),  as  explained  by  Dundurs  and  Santare  (5),  any  nonnegative,  locally 
integrable  function  f(x)  for  which 


f(x)dx  =  1  , 


IR 


and  such  a  function  f(x)  containing  (3,  where  (3  >  0, 


Mx)  =  (|*  - 


approaches  the  Dirac  delta  function,  <5(a;)  as  (3  approaches  zero: 


(16) 


(17) 


lim  fp(x)  =  6(x).  (18) 

This  technique  is  known  as  rnolliftcation,  where  the  function  fp(x)  is  a  mollifier  or  “averaging 
kernel”  and  any  function  that  is  convoluted  with  that  function  has  its  sharp  features  smoothed  (it 
is  mollified).  The  technique  is  used  here  to  provide  that  a  linear  differential  operation  on  the 
function  f(x)  produces  the  delta  function  in  accordance  with  equation  15,  and  as  such,  may  be 
thought  of  as  requiring  the  solution  for  the  concentrated  unit  of  inhomogeneity  to  be  a  distribution 
of  stresses  around  the  singularity  rather  than  an  infinite  stress  at  the  singularity.  Thus,  the 
singularity  is  smoothed  or  mollified. 

The  Dirac  delta  function  has  the  properties  of  being  zero  everywhere  except  at  the  point  of  some 
singularity,  say,  at  point  x  —  £,  where  the  Dirac  delta  function  suddenly  (conceptually) 
approaches  infinity.  Then  the  integral — over  any  domain,  say  R,  containing  that  point  of 
singularity,  x  =  £,  as  long  as  that  point  is  not  one  of  the  end  points  of  the  domain — of  any 
function,  such  as  equation  17,  with  these  properties  of  the  Dirac  delta  function  will  equal  one  as 
stated  by  equation  16. 

In  mathematics,  the  Green’s  function  is  defined  as  the  solution  to  an  inhomogeneous  differential 
equation  (e.g.,  the  boundary  value  problem  of  an  infinite  elastic  medium  containing  a  dislocation) 
that  is  subject  to  a  concentrated  unit  impulse  (e.g.,  the  delta  function  representing  the  concentrated 
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unit  of  inhomogeneity — the  dislocation)  at  a  prescribed  point  £  =  (£,  77)  e  £1  inside  the  domain 
(e.g.,  the  location  of  the  dislocation  in  an  infinite  elastic  medium).  The  Green’s  function  solution 
is  the  linear  superposition  (e.g.,  using  an  integral  in  a  continuum  to  accumulate  the  stress 
disturbances  due  to  the  dislocation  in  the  elastic  medium  concentrated  at  position  £)  to  sum  up 
Green’s  function  solutions  to  (homogeneous)  individual  unit  impulse  problems.  The  Green’s 
function  is  the  solution  to  the  homogeneous  portion  of  the  differential  equation  everywhere  x  ^  £. 
Then  the  solution  to  the  boundary  value  problem,  which  is  "forced"  and  so  a  non-homogeneous 
differential  equation  (e.g.,  not  to  confuse  with  the  term  "inhomogeneity"  used  above;  here  we 
mean  a  differential  equation  not  equal  to  zero,  p{x)^ |  +  q(x)^  +  r(x)u{x )  =  f(x),  where  f{x) 
is  the  forcing  function),  is  a  linear  superposition  of  the  Green’s  function  solutions  to  the  individual 
unit  impulse  problems  (11).  Then,  once  the  Green’s  function  is  known,  the  solution  to  the  forced 
boundary  value  problem  is  found  from  a  convolution  integral  of  the  forcing  function  (itself  a 
superposition  of  delta  impulses)  and  the  Green’s  function,  for  example,  in  two  dimensions 


u{x,y)=  /  /  G(x,  y\ £,  rj)f  (£,  rf)d^dr] .  (19) 

J  Jfi 

Further,  in  terms  of  dislocations,  taking  the  derivative  with  respect  to  x  of  the  Green’s  function  at 
the  point  of  concentrated  inhomogeneity  g(x ,  £),  i.e.,  (a  linear  differential  operation),  will 
result  in  a  function  with  the  properties  described  by  equations  16,  17,  and  18.  That  is,  the 
derivative  of  the  Green’s  function  must  have  a  jump  discontinuity  at  the  impulse  point  x  —  £.  So 
we  can  use  the  singularity  of  a  jump  (discontinuity)  in  displacement  seen  across  a  dislocation  as  a 
forcing  function  because  it  has  these  properties,  and  so  it  is  the  solution  to  the  boundary  value 
problem  of  the  otherwise  undisturbed  plane  containing  a  dislocation.  The  Green’s  function  then 
characterizes  how  the  system  responds  to  the  forcing  function  and  that,  in  terms  of  mechanics,  is 
the  stress  field  of  the  system.  See  Dundurs  and  Santare  (8)  for  further  discussion  on  using  a 
dislocation  as  the  singularity  that  provides  a  Green’s  function  that  gives  the  stress  field  solution  to 
the  problem. 

3.3  Dislocations  and  Cracks 

Further,  Dundurs  and  Santare  (5)  point  out  the  similarity  of  a  crack  to  a  series  of  dislocations. 
Then,  if  we  know  the  solution  for  a  single  dislocation,  this  opens  the  possibility  of  using  a 
superposition  of  that  solution  for  the  solution  of  a  pile-up  of  dislocations.  If  we  know  the 
solution  for  a  pile-up  of  dislocations,  which  is  similar  mathematically  to  a  crack,  that  is,  a  singular 
disturbance  in  an  elastic  body,  then  we  can  use  this  solution  as  the  Green’s  function  solution  to 
the  plane  elasticity  boundary  value  problem  of  interest.  Bilby  and  Eshelby  (12),  in  1968,  were 
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the  first  to  discuss  the  use  of  a  series  of  dislocations  as  the  solution  for  the  stress  field  of  a  crack  in 
an  elastic  medium.  The  following  discussion  is  based  on  the  discussion  by  Dundurs  and 
Santare  (5),  which  is,  in  turn,  based  on  the  discussion  by  Bilby  and  Eshelby  (72). 

Consider  an  infinite  linear  elastic  medium,  which  contains  a  continuous  line  of  dislocations  along 
the  x-axis,  —  a  <  x  <  a,  as  seen  in  figure  7.  If  we  were  to  magnify  our  view  of  the  pile-up  of 
dislocations,  it  might  look  something  like  figure  8.  In  figure  8,  we  see  that  pairs  of  each  of  three 
identical  edge  dislocations  adds  an  incremental  step  in  displacement  of  thickness  b  to  a  region 
(the  crack).  If  we  add  more  and  more  dislocation  pairs  and  make  the  dislocations  smaller  and 
smaller  in  thickness  b,  then  in  the  limit  we  approach  a  continuous  distribution  of  dislocations 
along  the  x-axis  that  may  be  represented  by  some  function  B(x). 


y 


Figure  7.  Infinite  elastic  body,  loaded  at  infinity  and 
containing  a  crack  of  length  2a. 


We  are  interested  in  knowing  the  stress  field  around  a  crack.  The  crack  lies  in  an  infinite  elastic 
medium  (figure  7)  and  we  know  the  size  of  the  crack  (  —a  <  x  <  a)  and  we  know  the  elastic 
material  properties:  Poisson’s  ratio  v,  shear  modulus  /x,  and  the  relationships  in  equation  20. 
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Figure  8.  Schematic  interpretation,  from  Dundurs  and  Santare  (8),  of  what 
a  crack  represented  by  dislocations  might  look  like. 


!3  —  4 v  plane  strain 
| plane  stress 

The  crack  is  represented  by  a  series  of  dislocation  pairs  having  a  distribution  that  is  given  by  a 
forcing  function  B(x,  y ).  Then,  to  determine  the  stress  field  around  the  crack  in  the  elastic 
medium,  we  simply  add  up  all  the  little  bits  of  stress  caused  by  all  the  little  bits  of  deformation 
caused  by  all  the  little  dislocations  described  by  B(x,  y )  with  the  point  of  singularity  of  interest 
(i.e.,  the  crack  tip  or  dislocation)  located  at  (£,  rj)  .  That  is, 


ayy(x,y)  =  J  B(Z,v)^<7w(.x,y,£,ri)d£dr} ,  (21) 

and  the  stress  components  in  the  other  directions  are  similar,  for  example,  with  crxx  or  axy 
replacing  cryy.  In  equation  21,  we  have  an  integral  equation  of  the  Green’s  function  solution, 

G  =  j;(Tyy(x,  y\  £,  rj),  for  a  single  edge  dislocation  located  at  point  (x,  y)  =  (£,  rj),  with  the  forcing 
function  B{x). 

So  if  we  consider  the  continuous  distribution  of  dislocations  to  lie  on  the  x-axis,  where  y  =  0  so 
that  the  dislocation  distribution  may  be  written  as  B(x),  and  consider  the  dislocation  distribution 
to  lie  between  —a<x<a  and  have  a  Burgers  vector  in  the  y  direction,  then  we  can  see  that  at 
some  location  between  x  and  x  +  dx  anywhere  along  this  line,  the  net  change  in  dislocations  is 
given  by  db  =  B{x)dx.  This  simplifies  equation  21  to 


tfyfav) 


2/i 

7 t(k  +  1) 


dj, 


—  a  <  x  <  a  . 


(22) 
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The  preceding  discussion  is  true  for  Mode  I  crack  opening  (i.e.,  tension)  and  the  stress 
components  in  the  other  directions  are  found  in  a  similar  manner.  We  could  also  find  the  stress 
components  for  the  other  crack  opening  modes,  e.g..  Mode  II  (i.e.,  anti-plane  shear)  using  screw 
dislocations,  and  similarly  Mode  III,  using  appropriate  orientations  in  each  case  for  the 
dislocation  distribution  and  corresponding  Burgers  vector. 

Similar  to  the  discussion  of  inclusions  above,  if  we  observe  that  figure  8  represents  a  pile-up  of 
tiny  steps  in  displacement,  then  if  we  insert  the  accumulation  of  those  jumps  in  displacement 
(i.e.,  the  crack)  into  an  undisturbed  elastic  medium,  we  would  cause  some  state  of  stress  in  the 
elastic  medium.  Thus,  for  our  collection  of  dislocations  along  the  x-axis,  where  —a  <  x  <  a  and 
y  —  0,  tractions  would  be  instigated  in  the  otherwise  undisturbed  elastic  medium.  Thus,  our  now 
disturbed  elastic  medium  would  experience  compressive  tractions  along  the  x-axis,  where 
—a  <  x  <  a  and  y  —  0  due  to  the  inserted  accumulation  of  displacements,  and  it  would 
experience  tensile  tractions  everywhere  else  along  the  x-axis,  that  is,  |x|  >  a. 

In  fracture  mechanics,  we  typically  impose  boundary  conditions  requiring  that  the  crack  faces  be 
traction  free  (13,  14).  But  we  just  realized  that  there  is  a  compressive  traction  in  the  region  of  the 
crack  due  to  the  inserted  continuous  distribution  of  dislocations,  and  we  found  that  equation  22 
describes  the  stress  in  the  region  due  to  this  insertion.  If  we  realize  that  we  stated  earlier  that  the 
infinite  elastic  medium  was  undisturbed  such  that  there  are  no  forces  applied  at  infinity,  then  we 
can  see  that  if  we  were  to  apply  some  stress  at  infinity,  say,  ayy,  we  could  counteract  the  stress  due 
to  the  crack,  ayy,  by  ensuring  that  ayy  +  ayy  =  0.  Then  if  we  were  to  remove  the  pile-up  of  tiny 
steps  in  displacement,  which  we  know  are  described  by  the  function  B(x),  then  that  same  space 
would  be  empty,  like  a  crack,  but  it  is  still  described  by  the  function  B(x)  and  the  stress  field 
around  it  is  still  described  by  ayy.  Thus,  we  have  described  the  crack  in  terms  of  the  tiny  steps  in 
displacement  caused  by  a  pile-up  of  dislocations. 

Typically,  we  know  or  can  measure  the  value  of  the  forces  applied  at  infinity,  oyy,  but  we  do  not 
know  the  function,  B(x,  y),  that  describes  the  continuous  distribution  of  dislocations  (a.k.a.,  the 
crack).  So  if  we  have  a  specimen  containing  a  straight  crack  of  size  2a  and  we  measure  the  forces 
applied  at  infinity,  then  we  let  the  coordinate  axes  lie  on  the  crack  such  that  the  crack  lies  along 
the  x-axis  with  the  center  at  y  =  0  and  —a  <  x  <  a,  as  in  figure  7,  then  we  can  find  the  unknown 
function  B(x ,  0)  using  equation  22  set  equal  to  the  negative  of  the  tractions  applied  at  infinity,  ayy. 
Then  if  we  consider  the  point  of  singularity,  the  crack  tip,  which  lies  at  (x,  y)  =  (£,  rf),  we  have 

o /„ = -< ’  -“<*<“■  <23) 
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but  we  need  another  constraint  to  solve  this  equation.  If  we  recognize  that  where  the  crack  closes 
there  is  no  step  in  displacement,  then  we  have  the  crack  closure  condition, 


B(£)dl £  =  0.  —  a  <  x  <  a 


The  solution  to  equations  23  and  24  is  well  known  and  is  given  by 


B(x)  =  —a, 


a(«  +  1)  x 

yv  2/x  y/a2  —  x2  ’ 


a  <  x  <  a  , 


(24) 


(25) 


B{x)  =  0  ,  \x\  >  a  ,  (26) 

Equations  25  and  26  simply  state  that  the  function  that  describes  the  continuous  distribution  of 
dislocations  is  equal  to  equation  25  on  the  crack,  but  it  is  equal  to  zero  away  from  the  crack  where 
there  is  no  displacement  due  to  dislocations. 

Now  we  know  the  function  B(x,  0)  that  will  result  in  zero  tractions  on  the  crack  surface.  That  is, 
cryy(x,  0)  =  0  on  —a  <  x  <  a.  Likewise,  due  to  the  tractions  at  infinity,  the  stresses  in  the 
specimen  everywhere  else  on  x,  along  y  —  0,  away  from  the  crack,  will  be  given  by 


ayy(x,0) 


(27) 


In  similar  fashion,  solutions  may  be  found  for  Mode  II  and  Mode  III  cracks  as  described  by 
Dundurs  and  Santare  (5).  Also,  now  we  have  the  solution  for  an  infinite  elastic  matrix  containing 
a  crack,  which  is  represented  by  dislocations.  The  solution  for  the  interaction  of  a  crack  with  an 
inclusion  then  may  be  found  by  using  the  solution  from  the  interaction  of  an  edge  dislocation  with 
an  inclusion  as  the  Green’s  function  in  the  problem  of  a  crack  interacting  with  an  inclusion  in  an 
infinite  elastic  medium. 


3.4  Green’s  Function  for  a  Dislocation 


In  an  infinite  elastic  body  with  no  tractions  applied  at  the  boundaries  (at  infinity),  for  a  single 
edge  dislocation  located  within  this  body,  the  elastic  stress  field  in  the  body  due  to  this  dislocation 
is  given  by  the  following  (Cartesian)  stress  components, 
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Vxx(x,y,t„v) 


2 yby  (x  -  0  [(s  -  O2  -(y-  rj)2\ 
tc(k  +  1)  r4 


(28) 


Vyy(x,y,£„y) 


2 yby  (x  -  Q  [(s  -  Q2  +  3 (y  -  rj)2] 
7 t(k  +  1)  r4 


(29) 


Vxy(x,y,t„y) 


2fiby  (y  -  r])[(x  -  £)2  -  (y  -  ri)2] 
7 t(k  +  1)  r4 


(30) 


where  r2  =  (x  —  £)2  +  (y  —  ;/)2,  by  is  the  Burgers  vector  oriented  in  the  y  direction,  //  is  shear 
modulus,  and  n  is  a  relationship  to  Poisson’s  ratio  as  described  in  equation  20  (5).  The  stresses 
may  be  found  for  any  arbitrary  location  in  the  Cartesian  plane,  say,  a  distance  of  r  from  the  origin, 
at  an  angle  of  6  from  the  x-axis,  using  a  simple  transformation  of  axes,  for  example  (15), 


alx  —  alx  cos2  0  +  aly  sin2  9  —  2 a].y  sin  9  cos  9  , 

(Jyy  =  oxxx  sin2  9  +  ( Tyy  cos2  9  +  2cr^  sin  9  cos  9  ,  (31) 

°xy  =  ialx  -  ^ )y )  sin  ^  COS  0  +  (cOS2  0  -  sill2  9)  . 

So  if  we  let  the  edge  dislocation  be  at  the  origin  of  the  coordinate  axes,  that  is,  (£,  rj)  =  (0, 0),  and 
then  if,  like  a  crack  along  the  x-axis  at  y  =  0,  we  look  closer  and  closer  at  the  dislocation  (which 
is  oriented  in  the  y-direction)  so  that  in  the  limit  the  dislocation  approaches  the  x-axis  as  y  — >  0, 
then  from  equation  29  we  have 
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<7yy(x,y,  °,0)  = 


•  •  ayy  oc 


2 fiby  (x  -  0)  [(x  -  0)2  +  3 (y  -  0)2] 
7r(/c  +  l)  ((x  -  0)2  +  (y  -  0)2)2 
/ A\  x(x2  +  3 y2) 

W  ( x2+y 2)2 


x 

(. x 2  +  y 2) 


2  y2 

(x2  +  y2) 
'l/x2\ 

Mx2) . 


2  y2  \ 

(x2  +  y 2) ) 


V  ( x 2  +  y2) 

1  l/x 

71  (i  +  (y/x)2) 


/A\ 

l/x 

V71"  J 

(1  +  yy/x)  ) 

(32) 


where  A  =  is  a  constant  dependent  on  material  properties  and  is  introduced  here  only  for 
convenience.  Taking  the  limit  of  ayy  as  y  0  we  see  that 


1  l/x 

lim  oc  —  - 77- 

y^ o  yy  it  (i  +  (y/x)  ) 

lim  ayv  (X  —  , 

y-> o  x 

so  the  stress  is  singular  (that  is,  a  singularity  or  a  point  force)  at  the  point  of  the  dislocation.  Also 

if  we  let  x  =  (3,  where  u  is  any  number,  then  take  the  limit  of  ayy  as  (3  — *  oo  we  see  that 


a yy{xi  y) 
&yy{X-)  y) 


x=0 


x=/3 


OC 


OC 


1  VI 

vr  (1  +  (y/P)2) 

1  P 

it  (1  +  P2y2) 


(33) 


so  that 
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0,  0 
(X)  ,  y  =  0 


(34) 


jim  «»)  = 


and 


1  P 

7r  (1  +  f32y 2) 


—  tan  1  f5y 

71 


oo 


y=- oo 


i 

—  7T  =  1 
7T 


(35) 


Notice  that  equations  34  and  35  satisfy  the  requirements  for  the  delta  function  given  by 
equations  12  and  13  (note  that  equation  13  is  the  same  requirement  as  given  by  equation  16).  So 
we  can  conclude  that  cryy  satisfies  the  requirements  for  the  delta  function. 

Also  notice  that  if  we  evaluate  the  limit  in  equation  34  as  (3  — *  0  we  see  that 


<7. 


yy 


oc  6{y) 


>o+ 


(36) 


which  is  the  same  as  concluded  by  equation  18.  Similar  arguments  may  be  made  for  the  other 
stress  components.  So  the  al3  satisfy  the  requirements  for  a  delta  function  and  so  may  be  used  as 
the  Green’s  function  in  integral  equations. 

Therefore,  we  now  know  that  the  Green’s  function  G(x,  y.  £,  rj)  as  in  equation  19  is  given  by  al3 
for  a  dislocation  interacting  with  a  boundary,  and  al3  satisfies  the  requirements  of  a  delta  function 
so  that  we  can  use  the  properties  of  equation  14  to  solve  equation  21  for  the  unknown  forcing 
function  B(x).  The  forcing  function,  call  it  B(x)  in  equation  21  or  call  it  /(£,  rj)  in  equation  19, 
or  whatever,  is  known  for  the  problem  of  the  stress  field  in  an  infinite  elastic  medium  caused  by 
the  presence  of  a  crack,  which  is  represented  as  a  pile-up  of  dislocations  as  discussed  previously. 
But  now  we  know  how  to  find  the  forcing  function  for  the  problem  of  the  stress  field  in  an  infinite 
elastic  medium  caused  by  the  presence  of  a  crack  near  an  inclusion.  All  we  need  are  the  Green’s 
functions  for  the  integral  equations  and  we  can  solve  the  integral  equations  for  the  forcing 
functions.  We  know  ai3  can  be  used  as  the  Green’s  functions,  so  we  use  equations  28,  29,  and  30. 
Using  these  stress  solutions  as  our  Green’s  functions,  we  can  find  the  unknown  forcing  functions. 
When  we  have  the  forcing  functions,  we  can  solve  for  the  stress  field  around  the  crack  and 
inclusion.  From  the  stresses,  we  can  determine  the  stress  intensity  factors  near  the  crack  and  the 
inclusion,  which  tell  us  how  the  crack  will  behave  and  are  important  for  predicting  fracture.  Also 
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once  we  know  the  stresses,  we  can  describe  the  strains  by  a  constitutive  model  such  as  Hooke’s 
law.  We  can  then  also  find  the  displacements  by  integrating  the  strains. 


4.  Solving  the  Crack-Inclusion  Interaction 


The  first  published  solution  to  the  crack-inclusion  interaction  problem  is  attributed  to  Tamate  in 
1968  (16).  Tamate  provided  the  solution  for  the  effect  of  a  circular  elastic  inclusion  on  a  crack. 
He  found  the  stress  field  around  a  crack  and  an  inclusion  mathematically  using  the  complex  stress 
potential  method;  a  famous  example  of  this  method  is  the  Westergaard  solution  for  the  stress  field 
around  a  crack  in  an  infinite  elastic  medium  (13).  Sanford  provides  a  review  of  solving 
mechanics  problems  using  complex  potentials  (14).  Tamate  followed  Muskhelishvili’s  method 
for  the  solution  of  complex  potentials  using  a  power  series  expansion  to  determine  a  set  of 
simultaneous  equations  to  determine  the  coefficients  in  the  power  series  (17).  In  his  example, 
Tamate  solved  20  equations  in  20  unknowns.  A  major  shortcoming  of  the  Tamate  solution  is  that 
the  solution  only  works  when  the  radius  of  the  inclusion  is  much  less  than  the  distance  from  the 
center  of  the  inclusion  to  the  crack,  and  then  only  when  that  distance  (center  inclusion  to  the 
crack)  was  much  less  than  that  distance  plus  the  crack  length,  that  is,  for  relatively  small 
inclusions  near  to  relatively  long  cracks.  This  shortcoming  was  pointed  out  and  addressed  by 
Atkinson  in  his  solution  for  the  crack-inclusion  problem  (18).  However,  the  Tamate  paper,  in 
addition  to  being  the  first  solution,  provided  a  critical  finding  for  the  fracture  mechanics  of  the 
crack-inclusion  interaction:  an  inclusion  more  rigid  than  the  surrounding  matrix  will  decrease 
the  stress  intensity  factor  in  the  neighborhood  of  the  crack  tip,  and  likewise  an  inclusion  less  rigid 
than  the  surrounding  matrix  will  increase  the  stress  intensity  factor  in  the  neighborhood  of  the 
crack  tip  (8,  16). 

Atkinson  in  1972  was  the  first  to  consider  the  crack  as  a  continuous  distribution  of  dislocations  to 
develop  a  singular  integral  equation.  He  then  solved  the  integral  equation  by  two  independent 
methods:  first,  he  reduced  the  integral  equation  to  a  Fredholm  equation  and  solved  the  Fredholm 
equation  numerically  by  a  method  Atkinson  published  in  a  separate  paper  (18).  This  is  the 
approach  taken  in  this  report.  The  second  method  used  an  expansion  of  Tschebyscheff 
polynomials  (18).  In  setting  up  the  integral  equations,  Atkinson  was  also  the  first  to  use  the 
Dundurs  and  Mura  solution  for  an  edge  dislocation  interacting  with  a  circular  inclusion  as  the 
Green’s  function  for  the  solution  of  a  crack  interacting  with  an  inclusion  (7,  18). 

The  Atkinson  solution  considered  only  cracks  and  dislocation  pile-ups  that  lie  on  the  x-axis 
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(i.e.,  where  y  =  0),  which  is  a  commonly  used  simplification  in  fracture  mechanics  (13,  14). 
However,  shortly  after,  in  1974,  Erdogan  et  al.  (19)  used  a  similar  approach,  that  is,  the  Dundurs 
and  Mura  solution  as  a  Green’s  function  in  a  set  of  singular  Fredholm  integral  equations,  to  solve 
the  more  general  crack-inclusion  interaction  for  an  arbitrarily  oriented  crack.  The  technique  used 
by  Erdogan  et  al.  has  since  become  the  most  widely  accepted  and  used  method  in  the  solution  of 
various  crack-inclusion  problems;  Dundurs  and  Santare  provide  a  much  more  comprehensive 
review  of  the  literature  in  this  field  (5). 

4.1  Forming  the  Integral  Equations 

An  integral  equation  includes  an  unknown  function  <p(t)  that  we  want  to  find  and  it  is  related 
through  an  integral  of  a  known  kernel  and  the  unknown  function.  If  both  limits  are  prescribed, 
we  call  the  equation  a  Fredholm  integral  equation  (if  one  is  variable,  it  is  called  a  Volterra  integral 
equation)  (20).  For  example,  a  Fredholm  integral  equation  of  the  first  kind  is 

f(x)=  f  K(x,t)cf>(t)dt  (37) 

J  a 

where  the  function  f(x)  and  the  integral  kernel  K(x,  t)  are  known  and  we  want  to  determine  the 
function  0(f).  Similarly,  a  Fredholm  integral  equation  of  the  second  kind  is 

<j>(x)  =  f(x)  +  A  f  K(x,t)4>(t)dt  (38) 

J  a 

where  A  is  some  known  constant  (usually  based  on  the  elastic  constants  of  our  materials)  and,  as 
before,  the  function  f(x)  and  the  integral  kernel  K(x,  t)  are  known  and  we  want  to  determine  the 
function  0(f),  but  now  we  also  have  a  function  <p(x).  The  function  (f>(x)  is  often  found  using  a 
Neumann  series  (20). 

We  know  from  the  earlier  discussions  (see  sections  3.2  and  3.4)  that  the  stress  field  solution  for  an 
edge  dislocation,  equations  28-30,  has  the  properties  of  the  delta  function  after  being  operated  on 
by  a  linear  differential  operator  and  so  may  be  used  as  a  Green’s  function  according  to 
equation  15.  We  know  that  a  crack  may  be  represented  mathematically  as  a  continuous 
distribution  of  dislocations,  and  we  know  that  we  can  determine  the  stress  field  of  a  crack 
interacting  with  an  inclusion  if  we  have  the  Green’s  function  (see  section  3.3).  So  we  can  use  the 
stress  field  solution  as  the  Green’s  function  in  the  integral  equations  that  describe  the  stress  field 
solution  for  a  crack  interacting  with  an  inclusion.  So  starting  simply,  with  equation  37,  we  let  the 


21 


kernel  K(x,  t )  equal  the  Green’s  function  for  an  edge  dislocation  given  by  equations  28-30. 

Presumably,  we  also  know  the  size  of  the  inclusion  and  the  size  of  the  crack,  so  as  we  build  our 
integral  equations,  we  set  the  limits  of  integration  to  the  limits  of  inclusion  size  and  crack  size. 
Presumably  also,  we  know  the  stress  field  in  our  infinite  elastic  medium.  We  may  find  that  we 
have  other  functions  in  our  integral  equations  and  if  we  bring  constants  out  in  front  of  the  integral, 
for  example,  as  A,  we  may  find  that  our  integral  equation  looks  more  like  equation  38.  Finding 
the  solution  to  our  integral  equations  will  involve  numerical  methods. 

4.2  Inclusion  and  Crack  Geometry 

Solutions  have  been  found  for  inclusions  that  are  flat,  circular,  or  elliptical,  and  cracks  may  be 
straight  or  arc  and  cracks  may  be  present  in  the  matrix  or  the  inclusion  or  both.  Although  a 
general  solution  is  provided  by  Erdogan  et  al.  for  an  arc  crack  (19),  they  necessarily  simplify  to  a 
straight  crack  to  avoid  complicating  coordinate  system  transformations.  For  simplicity,  here  we 
consider  only  a  single  straight  crack  in  an  infinite  elastic  medium.  Further,  we  assume  that  the 
crack  and  inclusion  are  about  the  same  order  of  magnitude.  We  work  in  two  dimensions  and 
assume  that  inclusions  are  “sparsley”  located  since,  in  the  words  of  Erdogan  et  al.,  “The  exact 
elasticity  treatment  of  the  three-dimensional  problem  with  a  regular  or  random  array  of  elastic 
inclusions  embedded  into  an  elastic  matrix  containing  an  internal  crack  appears  to  be  hopelessly 
complicated.  (79)”  Other  simplifying  assumptions  that  we  make  include  the  use  of  plane  strain  or 
generalized  plane  stress,  and  we  require  a  perfect  bond  between  the  inclusion  and  the  matrix. 
Finally,  we  consider  a  flat  inclusion.  The  geometry  we  have  just  described  for  the  problem  of 
determining  the  stress  field  near  a  crack  and  inclusion  in  an  infinite  elastic  medium  is  shown  in 
figure  9.  Here,  we  have  chosen  to  use  the  notation  provided  by  Erdogan  and  Wei  (21).  Note  that 
our  inclusion  is  specified  by  hQ(x),  but  that  since  we  are  assuming  a  flat  inclusion  we  require  ha  to 
be  a  constant  thickness  for  the  inclusion  and  ha  must  be  small  relative  to  the  length  of  the 
inclusion  2a\.  However,  we  could  specify  a  function  to  describe  the  inclusion  as  a  circle  or  an 
ellipse;  for  example,  if  we  assume  the  size  2a!  is  fixed  for  the  inclusion  (so  that  a  and  b  are 
constant  limits  for  our  integral  equations),  then  we  could  specify  a  function,  say 
hQ(x)  =  bQ  Vl  —  x2,  to  provide  an  elliptical  inclusion.  Finally,  as  an  “infinite”  elastic  medium, 
we  assume  that  the  boundaries  of  the  medium  are  sufficiently  far  away  from  the  region  of  the 
crack-inclusion  interaction  that  the  boundaries  do  not  experience  the  stress  influence  from  the 
crack  or  inclusion. 

In  figure  9  we  have  a  flat  inclusion  of  size  2 ax  with  limits  a  and  b  and  constant  height  hQ  and 
which  lies  on  the  x\,  y\  coordinate  axes.  We  also  have  a  straight  crack  of  size  2a2  with  limits  c 
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Figure  9.  Geometry  for  the  crack-inclusion  problem,  adapted  from 
Erdogan  and  Wei  (21) . 


and  d  that  lies  on  the  x2, 2/2  coordinate  axes.  The  x2,  y-i  coordinate  axes  are  oriented  relative  to 
the  x  1 ,  yi  coordinate  axes  by  a  rotation  of  angle  0.  Far  from  the  crack-inclusion  interaction 
region,  at  the  boundaries  of  the  infinite  elastic  medium,  we  have  applied  tractions  producing 
Cartesian  stresses  cr°?  as  shown  in  figure  9. 

4.3  Displacement  Derivatives 

We  have  defined  the  geometry  and  the  boundary  conditions  of  our  problem.  We  can  assume  that 
we  know  the  elastic  constants  /j  and  v  (and  from  v  we  have  k,  given  by  equation  20)  for  the  matrix 
and  inclusion  (or  void,  where  fi  =  0).  We  have  identified  the  Green’s  function  for  our  integral 
equations.  Now  we  need  to  identify  the  forcing  function  to  complete  our  integral  equations. 

Using  similar  logic  to  the  earlier  discussions,  if  we  have  an  undisturbed  matrix  and  we  insert  an 
inclusion  or  we  insert  a  crack,  we  would  be  inserting  volume  that  would  have  a  displacing  effect 
on  the  surrounding  matrix.  So  we  use  the  jumps  in  displacement  across  the  inclusion  and  across 
the  crack  as  our  concentrated  units  of  inhomogeneity.  For  rq  and  vt,  which  are,  respectively,  the 
Xi  and  yi  components  of  displacement  for  i  —  1,2,  the  coordinate  systems  shown  in  figure  9. 

Then  we  take  the  derivative  of  these  jumps  in  displacement  and  define  our  unknown  functions: 
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(39) 


9i(xi)  =  [t»i(xi,  0+)  -  vi(x!,  0  )]  ,  (a  <  X!  <  b)  , 

Q 

h\(x\)  =  —  [mi(xi,0+)  -mi(xi,0-)]  ,  (. a  <  xi  <  b)  ,  (40) 

d 

g2(x2)  =  7^-  [v2(x2,  0+)  -  v2(x2,  0')]  ,  (c  <  Xx  <  d)  ,  (41) 

Q 

h2(x2)  =  —  [u2(x 2,  0+)  -  M2(a;2,  0  )]  ,  (c  <  Xi  <  d)  .  (42) 

We  want  to  know  what  the  stress  field  in  the  matrix  near  the  crack  and  inclusion  is.  Referring  to 
figure  9,  we  now  have  an  unknown  function  g\ ,  which  is  the  derivative  of  the  discontinuity  in  the 
displacement  of  the  matrix  in  the  y\  direction  across  the  inclusion  where  ij\  —  0  and  (a  <  X\  <b). 
Likewise,  we  have  an  unknown  function  h\,  which  is  the  derivative  of  the  discontinuity  in  the 
displacement  of  the  matrix  in  the  x\  direction  along  the  inclusion  where  y\  —  0  and  (a  <  x\  <  li). 
We  also  have  an  unknown  function  g2,  which  is  the  derivative  of  the  discontinuity  in  the 
displacement  of  the  matrix  in  the  y2  direction  across  the  crack  where  y2  =  0  and  (c  <  x2  <  d). 
Likewise,  we  have  an  unknown  function  h2,  which  is  the  derivative  of  the  discontinuity  in  the 
displacement  of  the  matrix  in  the  x2  direction  along  the  crack  where  y2  =  0  and  (c  <  x2  <  d). 

4.4  Stresses  in  the  Inclusion 

Considering  boundary  conditions  for  a  moment,  we  know  crack  faces  are  traction  free  so  we 
realize  that  the  stresses  must  be  zero  inside  the  crack,  but  the  displacement  jump  cannot  be  zero. 
Likewise,  outside  the  crack,  the  stresses  are  not  zero,  but  the  displacement  jump  is  zero.  On  the 
other  hand,  the  inclusion  is  not  free,  rather  we  assumed  it  is  perfectly  bonded  to  the  matrix,  and  let 
us  suppose  that  the  inclusion  has  some  known  elastic  constants  (i.e.,  it  is  not  void),  say,  E0,  ua,  /j0. 
Then  we  must  realize  that  the  stresses  both  inside  and  outside  of  the  inclusion  will  not  be  zero, 
and  likewise  the  jumps  in  displacement  will  not  be  zero  either  inside  or  outside  of  the  inclusion. 

Let  us  assume  that  the  variation  of  stress  along  the  length  of  the  inclusion  is  negligible  and 
likewise  the  strain  along  the  length  of  the  inclusion  is  negligible.  These  are  valid  assumptions 
because  the  thickness  of  the  inclusion  is  small  relative  to  its  length  and  because  the  length  of  the 
inclusion  will  only  change  by  a  negligible  amount  under  stress.  Let  us  also  consider  plane  strain 
conditions  here  so  that  for  the  elastic  constants  we  have 
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E0  =  2/x0(l  +  u0)  ,  k0  =  3  -  Av0  . 


(43) 


We  know  from  elementary  elasticity  theory  (15)  that  strain  is  equal  to  the  derivative  of 
displacement  or  even  more  basically  that  strain  is  equal  to  the  change  in  length  divided  by  the 
original  length.  That  is,  in  two  dimensions,  the  strains  in  the  x  and  y  directions  are 


du 

dv 

du 

dv 

dx  ’ 

£yy=dy ’ 

&xy  W  f 

dy 

dx 

(44) 


or  written  another  way, 


£k  = 


A Lk 

Lk 


(45) 


where  k  =  x,  y  are  the  x  and  y  directions. 

Then  let  us  observe  that  the  change  in  the  length  of  the  inclusion,  A  Lk  is  the  same  as  the  jump  in 
displacement  in  a  given  direction  (i.e.,  k  —  x,  y),  and  also  observe  that  the  original  length  in  the  y 
direction  is  h0(x).  Recall  that  we  assumed  that  the  change  in  length,  du,  is  negligible  along 
a  <  X\  <  b  so  that  we  know  from  equation  44  that  the  strain  in  the  inclusion  in  the  x  direction  is 
£xx  =  0-  Likewise,  since  du  =  0,  the  component  of  the  shear  strain,  exy  from  equation  44  that  is 
dependent  on  du,  is  also  zero. 

Therefore,  we  have  for  the  strains  in  the  inclusion 


Cc(*i)  = 


[t>i(a;i,0+)  -  ni(xi,0  )' 
h0{x  i) 


n  inc(  \  ^  [zZj  (q^i ,  0+)  —  Ui(x\,  0  )' 

XV  1  lJ  h0(Xl) 

Also,  from  equations  39  and  40,  we  have  the  partial  derivatives 

d  [i>i(xi,  0+)  -  vi(Xl,  0-)]  =  gi(Xi)dxi  , 


(46) 


(47) 


(48) 
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d  u\{xi,  0+)  —  u\(x\,  0  )  =  hi(x i)dxi  , 


(49) 


and  if  we  integrate  these  along  the  inclusion  we  have  an  expression  for  the  jumps  in  displacement: 


^1(^1,  0+)  —  Vi(xi,  0  ) 


J  gi(x1)dx1  , 


(50) 


Ml(xi,0+)  -Mi(xi,0  ) 


J  h\{x\)dx\  . 


(51) 


Hooke’s  Law  tells  us  that  stress  is  related  to  strain  by  way  of  the  elastic  constants  for  the 
material  (75),  so  we  can  write 


"yy 


(xi)  = 


Vn 


E0(  1 


2  7/  2 

o  me 

'  ayy 


v. 


(ah)  , 


(52) 


O 


£yy{xi)  =  {xi)  ’ 


(53) 


and  if  we  solve  for  the  stresses, 


=  •  (54) 

=  2^47(10.  (55) 

Then  if  we  plug  equations  50  and  5 1  into  equations  46  and  47  and  then  plug  the  result  into 
equations  54  and  55  and  finally  convert  the  elastic  constants  using  equation  43  (for  plane  strain), 
we  find  an  expression  for  the  stresses  in  the  inclusion  as  a  function  of  known  elastic  constants  and 
our  unknown  displacement  derivatives  defined  in  equations  39  and  40: 

/  91  (t)dt ,  (56) 

fi'o  -L  1 )  J  a 
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°Tyc(x i)  =  j^r y  fa  hi(t)dt  •  (57) 

4.5  Stresses  in  the  Matrix 

We  have  found  expressions  for  the  (in  plane)  stresses  in  the  inclusion  as  a  function  of  known 
elastic  constants  and  unknown  displacement  derivatives.  Next,  we  need  to  find  expressions  for 
the  stress  disturbance  in  the  matrix  due  to  the  crack  and  due  to  the  inclusion.  We  now  call  upon 
the  dislocation  solution  from  Dundurs  and  Mura  (7)  (or  similar).  Earlier  discussion  demonstrated 
that  we  can  use  the  dislocation  solution  as  the  kernel  in  a  set  of  integral  equations  to  determine  the 
stresses  in  the  plane  due  to  the  inclusion  and  the  crack.  Then  with  our  unknown  forcing  functions 
defined  earlier,  the  displacement  derivatives  of  equations  39-42,  we  have  the  components  of  a 
Fredholm  integral  equation  of  the  first  kind. 

The  in  plane  stresses  in  the  matrix  are  caused  by  the  presence  of  the  inclusion,  that  is,  caused  by 
the  jumps  in  displacement  due  to  the  presence  of  the  inclusion,  that  is,  caused  by  the  forcing 
functions  defined  in  equations  39  and  40,  the  displacement  derivatives.  These  equations  are  the 
4>(t)  functions  in  our  integral  equations  37.  The  dislocation  solution  given  by  equations  28-30 
characterizes  the  response  of  the  matrix  to  the  forcing  functions  from  the  inclusion.  So  the 
stresses  at  a  point  (aq,  y\ ),  referred  to  the  aq,  y\  coordinate  system,  due  to  the  inclusion  of  size 
a  <  x\  <b,  may  be  expressed  by 

all(xi,yi)=  [Gxx(xi,y!,t)gi(t)  +  Hxx(x!,yi,t)hi(t)]dt ,  (58) 

J  a 

<?w{xi,yi)=  [Gyy(x1,y1,t)g1(t)  +  Hvy(x1,y1,t)h1(t)]dt ,  (59) 

J  a 

<Tll(xi,yi)=  [GXy(xi,  2/1,  t)si(t)  +  Hxyfruyutfhiit)]  dt .  (60) 

J  a 

The  stresses  at  a  point  (x2, 1/2)7  referred  to  the  x 2,1)2  coordinate  system,  due  to  the  crack  of  size 
c  <  x2  <  d,  may  be  expressed  by 

fd 

<xll(x2,y2)=  [Gxx(x2,y2,t)g2(t)  +  Hxx(x2,y2,t)h2(t)]  dt ,  (61) 
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fd 

°H(x2,y2)=  [Gyy(x2,y2,t)g2(t)  +  Hyy(x2,  y2,t)h2(t)\  dt ,  (62) 


9  fd 

axy(x2i  y-i)  =  J  [Gxy{x2,y2,  t)g2(t )  +  HXy(x 2,  y2,  7)^00]  dt .  (63) 

In  the  tensor  notation  used  here  for  stress,  cr^- ,  we  have  i,  j  —  x,y  representing  the  Cartesian 
stress  components  in  the  usual  sense,  and  k,  l  —  1,  2,  where  k  =  1,  2  indicates  either  stress  caused 
by  the  inclusion  k  =  1  or  stress  caused  by  the  crack  k  =  2  and  (  =  1,2  indicates  stresses  referred 
to  either  coordinate  system,  xi,  y\  for  /  =  1  or  x2. 1)2  for  l  =  2.  Also  note  that  in  earlier  notation 
we  used  the  point  (x,  y)  =  (£,  rf)  in  our  integral  equations,  but  here  we  use  the  notation  provided 
by  Erdogan  and  Wei  (27)  so  that  the  point  (x,  y)  =  (£,  rf)  =  (t,  0),  where  since  we  assumed  a  flat 
inclusion,  we  have  y\  =  0  and  on  a  crack  we  have  y2  =  0. 

In  equations  58-63,  the  Green’s  functions  given  by  the  dislocation  solution  provided  in 
equations  28-30  (for  a  dislocation  oriented  in  the  y  direction)  must  be  rotated  for  edge 
dislocations  in  the  appropriate  x  or  y  direction,  that  is,  Gij  is  for  an  edge  dislocation  oriented  in 
the  y  direction  and  Hij  is  for  an  edge  dislocation  oriented  in  the  x  direction.  Note  that  x,  y  refer 
to  the  appropriate  coordinate  system,  either  x\,yi  or  x2,  y2,  and  /j  and  n  are  elastic  constants  of 
the  matrix  where  k  is  given  by  equation  20.  Then  we  have  the  following  Green’s  functions: 


Gxx(x,y,t,  0) 


2/i 


(1  +  k) 


(t  —  x)[(t  —  x)2  —  y2 
[(f  —  x)2  +  y2]2 


(64) 


Gyy[x,  y ,  t  ,  0) 


2/i 


L(i  +  /c) 


(t  —  x)[(t  —  x)2  +  3  y2 
[{t  —  x)2  +  y2]2 


(65) 


Gxy(x,y,t,  0)  = 


2/i 


[(!  +  «) 


(y)[-(t  -  x)2  +  y2 
[(f  —  x)2  +  y2]2  ’ 


(66) 


Hxx(x,y,t,  0)  = 


2/i 


{y)[y2  +  3(f  —  x) 


L(l  +  /c)J  [(t  —  x)2  +  y2]2 


(67) 


Hyy(x,y,t,  0)  = 


2/i 


(!  +  «) 


(: y)[y 2  -  (<  -  ^)2] 

[(f  -  x)2  +  y2}2  ’ 


(68) 
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Hxy(x,y,t,  0) 


(69) 


2  fj,  (■ t  —  x )  [— y2  +  (t  —  x )2] 

(1  +  k)  [(£  —  x)2  +  y2]2 

4.6  Coupling  Stresses  Between  Crack  and  Inclusion 

In  addition  to  having  an  influence  on  the  stress  field  in  the  matrix  referred  to  the  respective 
coordinate  systems  on  which  the  inclusion  or  crack  lies,  the  inclusion  influences  the  stress  in  the 
matrix  around  the  crack  and  the  crack  influences  the  stress  in  the  matrix  around  the  inclusion. 
These  are  coupling  stresses  between  the  inclusion  and  the  crack.  The  coupling  stresses  are  found 
using  a  transformation  of  coordinate  axes  as  in  equation  31,  where  the  angle  9  relates  the  two 
coordinate  axes  as  in  figure  9. 

Then  using  the  notation  described  above  for  afj,  the  coupling  stresses  in  the  matrix  may  be  found 
from  the  previously  defined  stresses  using  the  following  transformations.  Here  the  stresses  from 
equations  58-60,  which  are  calculated  at  point  (xi,yi)  above,  are  calculated  at  the  point 
(xi,yi)  =  (x2  cos  9,  x2  sin  9)  in  the  coupling  stresses  below,  and  likewise  the  stresses  from 
equations  61-63,  which  are  calculated  at  point  (x2, y2)  above,  are  calculated  at  the  point 
(x2,  IJ2)  =  (x1  cos  9,  —Xi  sin  9)  in  the  coupling  stresses  below. 

Recall  that  we  assumed  the  variation  in  stress  in  the  x  direction  is  negligible.  The  stress 
transformations  for  coupling  stresses  are  the  following: 

<?ll(x  1,  0)  =  cr22(x2,  y2)  sin2  9  +  cr^(x2,  y2)  cos2  9  +  2a2x2y(x2,  y2)  sin  9  cos  9  ,  (70) 


^5(^1,  °)  =  {^ll(x2,y2)  -  a2y2y(x2,y2))  sm9  cos9  +  a2x2y(x2,y2)  (cos2  9  -  sin2  9)  ,  (71) 

(JH(x2,  0)  =  alKx!,  yi)  sin2  9 +  aH(x1,yl)  cos2  9 +  2all(xi,yl)  sin  9  cos  9  ,  (72) 


2,0)  =  (crlKx  1, 2/1)  -  a^xx,  yi))  sin 9  cos 9  +  a^Xi,  (cos2  9  -  sin2  9)  .  (73) 
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Likewise,  the  Green’s  function  transformations  for  coupling  stresses  are  the  following. 

Plug  the  expressions  from  equations  58-63  into  equations  70-73,  ignore  for  now  the  integral 
signs  and  the  displacement  derivatives  (they  will  return  later),  and  collect  and  separate  G  and  H 
terms  to  find  the  Green’s  function  transformations  for  coupling  stresses  as  follows: 

Gfy(x i,  t)  =  Gxx(x2,  y2)  sin2  9  +  Gyy(x2,  y2)  cos2  9  +  2Gxy(x2,  y2)  sin  9  cos  9  ,  (74) 

Gl2y(xi,t)  =  (■ Gxx(x2,y2 )  -  Gyy(x2,y2))sm9cos9  +  Gxy(x2,y2)  (cos2  9  -  sin2  9)  ,  (75) 

Hyy(x  1,  t)  =  Hxx(x 2,  y2)  sin2  9  +  Hyy(x2,  y2)  cos2  9  +  2Hxy{x2,  y2)  sin  9  cos  9  ,  (76) 

t)  =  (Hxx(x 2,  y2)  -  Hyy(x2,  y2))  sin  9  cos  9  +  Hxy(x2 ,  y2)  (cos2  9  -  sin2  6>)  ,  (77) 

where  (x2,  y2)  =  (xi  cos  9,  —x\  sin  9),  and 

Gfy(x2,t)  =  Gxx (xi,  yi )  sin2  9  +  Gyy(xi,yi)  cos2  9  -  2Gxy(x1,y1)  sin9  cos  9  ,  (78) 

G2xy{x2,t)  =  (Gxa;(xi,yi)  —  Gra(xi,yi))  sin  0  cos  0  +  Gxy(xi,yi)  (cos2  6>  -  sin2  9)  ,  (79) 

H‘yy(x2lt)  =  Hxx(x1,y1)sin29  +  Hyy(x1,y1)cos29  -  2Hxy(x1,y1)sm9 cos9  ,  (80) 

H2l(x2,t )  =  (ifxx(xi,yi)  —  Hyy(xi,  yi))  sin  0  cos  0  +  Hxy(xi,  yi)  (cos2  9  -  sin2  61)  ,  (81) 
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where  (xi,yi)  =  (x2  cos9,  x2  sin#). 

Then  the  normal  and  shear  components  of  the  coupling  stresses  on  the  y\  —  0  plane  of  the 
inclusion,  which  are  due  to  the  influence  of  the  displacement  derivatives  g2(x2)  and  h2(x2)  for  the 
crack  are 


r d 

i,  0)  =  j I  [Gfy(x i,  t)g2(t)  +  H™(x i,  t)h2[t)\  dt ,  (82) 

i,  0)  =  j  [G™(x i,  t)g2(t)  +  H™{x i,  t)h2(t)]  dt ,  (83) 

and  the  normal  and  shear  components  of  the  coupling  stresses  on  the  y2  =  0  plane  of  the  crack, 
which  are  due  to  the  influence  of  the  displacement  derivatives  gi(x\)  and  h ,  (x  | )  for  the  inclusion 
are 


Vyl(x 2,0)  =  j  [G2yl(x2,t)g1(t)  +  H^(x2,t)h1{t)]  dt ,  (84) 

J  a 

°2xy(x 2,0)  =  f  [Gly(x2,t)gi(t)  +  H^y(x2,  t)hi(t)]  dt .  (85) 

J  a 

4.7  Boundary  Conditions 

Let  the  infinite  elastic  medium  be  uniformly  loaded  far  away  from  the  region  of  the  crack  and 
inclusion,  as  shown  in  figure  9.  If  the  xi,  i)\  coordinate  axes  are  aligned  with  the  direction  of  the 
a™x  and  stress  components  as  shown  in  figure  9,  and  recalling  again  that  stress  difference  in 
the  x  direction  was  assumed  negligible,  then  we  have  the  following  stresses  along  the  x\  axis  in 
the  inclusion: 


=  a7y  ’  (86) 

a1™(x1,0)  =  a™.  (87) 
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For  coordinate  axes  that  are  not  aligned  with  the  directions  of  the  uniform  loading  at  infinity,  such 
as  the  X2 , 1/2  coordinate  axes,  a  coordinate  transformation  as  performed  earlier  and  as  is  common 
in  elementary  elasticity  (15)  is  needed.  For  an  angle  6  between  the  direction  of  the  applied  loads 
at  infinity  and  the  coordinate  axes,  as  shown  in  figure  9,  a  transformation  is  performed  according 
to  equation  31.  Then  we  have  the  following  stresses  along  the  x2  axis  in  the  crack: 


(x2,  0)  =  sin2  6  +  a™  cos2  6  -  2a™  sin  d  cos  9  ,  (88) 

a2™  (x2,  0)  =  (a™  -  a™)  sin  6  cos  9  +  a™  (cos2  d  -  sin2  6)  .  (89) 

We  now  have  expressions  for  all  of  the  stresses  in  the  inclusion  and  in  the  matrix  around  the 
inclusion  and  the  crack  including  (1)  the  stresses  in  the  inclusion  caused  by  the  jump  in 
displacement  from  the  presence  of  the  inclusion,  which  pushes  on  the  matrix  and  causes  the 
matrix  to  push  back  on  the  inclusion  inducing  stress  within  the  inclusion;  (2)  the  stresses  in  the 
matrix  due  to  the  jumps  in  displacement  caused  by  the  presence  of  the  inclusion  and  the  crack;  (3) 
the  stress  caused  in  the  crack  by  the  jump  in  displacement  of  the  inclusion,  which  pushes  on  the 
matrix  and  causes  the  matrix  to  push  on  the  nearby  crack,  and  the  stress  caused  in  the  inclusion  by 
the  jump  in  displacement  of  the  crack,  which  pushes  on  the  matrix  and  causes  the  matrix  to  push 
on  the  nearby  inclusion;  and  (4)  the  stresses  caused  in  the  crack  and  in  the  inclusion  by  the 
tractions  applied  at  the  far  boundaries  of  the  infinite  elastic  matrix. 

If  we  recognize  that  for  equilibrium,  since  the  inclusion  is  assumed  to  be  perfectly  bonded  to  the 
surrounding  matrix,  the  stresses  in  the  inclusion  must  equal  the  stresses  in  the  surrounding  matrix, 
then  we  can  sum  the  stresses  in  the  matrix  around  the  inclusion  and  set  them  equal  to  the  stresses 
in  the  inclusion.  Then  we  have  the  following  boundary  conditions  in  the  inclusion  where  y\  —  0 
and  a  <  xi  <  b: 


<0^1, 0)  +  crUix^O)  +  al™(x  1,  0)  =  a^c(xi)  ,  (90) 

°)  +  (TS(aT>°)  +  (TS°(xi’°)  =  aZc(x  1)  •  (91) 

We  know  that  crack  faces  are  traction  free  (14),  so  inside  the  crack  there  must  be  zero  stress,  and 
so  if  we  sum  the  stresses  in  and  around  the  crack,  they  must  cancel  to  zero  so  that  there  is  no 


32 


stress  inside  the  crack.  Then  we  have  the  following  boundary  conditions  in  the  crack  where 

1/2  =  0  and  c  <  X2  <  d. 


<?fy(x 2,  0)  +  i Jyl(x2 ,  0)  +  a^(x2,  0)  =  0  ,  (92) 

2,  0)  +  afy(x i,  0)  +  al™(x2,  0)  =  0  .  (93) 

We  have  integral  expressions  for  each  of  the  components  of  equations  90-93  in  terms  of  known 
elastic  constants,  known  Green’s  functions,  and  unknown  displacement  derivatives.  Now  we  can 
construct  a  set  of  integral  equations  and  solve  for  the  unknown  functions. 

4.8  Integral  Equations  for  the  Crack-Inclusion  Interaction 

By  plugging  the  stress  expressions  given  by  equations  56  and  57,  59-63,  82-85,  and  86-89  into 
the  boundary  conditions  given  by  equations  90-93,  we  obtain  the  following  set  of  integral 
equations  for  the  crack-inclusion  problem. 

Along  the  inclusion,  where  ij\  —  0  and  a  <  X\  <  b,  we  have 


[Gyy(x1,t)g1(t)  +  Hyy(xi,t)hi(t)]dt  +  j  [Gjj (x  1 , t ) g2 (t )  +  H™(x  1,t)h2(t)]  dt  +  a\ 
fXl  ^0  +  1  I-1 


oo 

yy 


K0-  1  ho(x  i) 


9i(t)dt , 


(94) 


[Gxy(x1,t)g1(t)  +  HXy(x!,t)hi(t)]dt  +  /  [G™(xi,t)g2(t)  +  H™(xi,  t)h2(t)]  dt  + 


(7. 


x  y 


t*x  i 


ha{x  i) 


hi(t)dt . 


(95) 


Along  the  crack,  where  y2  =  0  and  c  <  x2  <  d, 
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(96) 


[Gyy(x2,t)g2(t)  +  Hyy(x2,t)h2(t)\  dt  +  [G 1 1 (x2 ,  t ) g  1  ( t )  +  Hyy(x2,t)h1(t)]  dt 

J  a 

+  cr“  sin2  9  +  cr“  cos2  0  —  2 a™  sin  6*  cos  9  =  0  , 


[Gxy{x2,t)g2(t)  +  Hxy(x2,t)h2(t)\dt+  /  [G2xy(x2lt)gi(t)  +  H^(x2,t)hi(t)]  dt 

J  a 

+  (cr~  -  cr~ )  sin  0  cos  9  +  cr“  (cos2  0  -  sin2  0)  =  0  . 


(97) 


These  equations  may  be  simplified  further  if  we  plug  the  equations  64-69  into  the  integral 
equations  94-97  where  appropriate  and  evaluate  the  relations  for  equations  64-69  at 
(xi,yi  =  xi,  0,  t,  0)  in  integrals  evaluated  along  the  inclusion  a  <  x\  <  b  and  evaluate  the 
relations  for  equations  64-69  at  (x2,  y2  =  x2, 0,  t,  0)  in  integrals  evaluated  along  the  crack 

c  <  x2  <  d. 

Then  from  equations  65  and  68  we  get 


Gyy{x,y,t,  0) 

Gyy^X i,  0,  t,  0) 


Gyy{xi,t ) 


2/i 


(t  —  x)[(t  —  x)2  +  3  y2 


(1  +  k)J 

[(t  —  x)2  +  y2}2 

2y 

( t  ~  xi)[(t  -  xi)2  +  3(0)2] 

(!  +  «). 

[(t~x  i)2 

+  (o)2]2 

2  /i 

(t  -  2fi)3]  _ 

2/i 

1 

(1  +  k). 

(t-x  i)4 

.(!  +  «). 

t  —  X  i 

(98) 


Hyy{x,y,t,  0) 

Hyy(X\^  0,  t,  0) 
Hyy(Xl,  0,  t,  0) 


2/i 

(l+«) 

2  fi 

(l  +  «) 


(y)[y2  -  (t  —  x)2] 

[{t  —  x)2  +  y2]2 

(0)  [o2  —  (t  —  ar)2] 

[(t  -  Xi)2  +  (0)2]2 


(99) 


and  plugging  these  into  equation  94  and  simplifying,  we  get 
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2/4  1  fb  9i(t)  dt  (1q(k0  +  1)  rxi  giit)  dt 

(1  +  K)irja  t-Xx  K0~  1  Ja  /l0(xi) 

rd  rd  ^  UU') 

+  J  G™(x1,t)g2(t)dt  +  J  H™(x1,t)h2(t)dt  =  -cr™  . 

Equation  100  may  be  further  expanded  by  plugging  in  equations  74  and  76,  which  are,  in  turn, 
expanded  by  plugging  in  equations  64-69.  Then  equation  101  results. 

This  expression  has  become  quite  cumbersome  so  we  do  not  expand  the  remaining  integral 
equations  completely  in  this  fashion;  this  expansion  is  conducted  here  only  to  show  the  complete 
integral  equation  for  the  normal  stresses  in  the  inclusion.  In  similar  fashion,  we  may  expand  and 
simplify  equations  95-97.  Thus  we  have  the  integral  equations  for  the  crack- inclusion  interaction 
problem.  Observe  that,  once  the  geometry,  as  in  figure  10,  and  material  properties  of  the  problem 
are  fully  defined,  the  integral  equations  are  composed  entirely  of  known  quantities  save  the  four 
unknown  displacement  derivatives,  equations  39-42.  So  we  have  four  integral  equations  in  four 
unknowns.  We  can  solve  for  the  four  unknowns  and  then  completely  describe  the  stress  field  of 
the  crack-inclusion  interaction.  However,  it  is  obvious  that  solving  these  four  integral  equations 
in  four  unknowns  is  only  practical  using  numerical  methods. 


2/i  1 

(1  +  k)  tt 


9l(t)  dt  _  do(^o  +  1) 
t  —  X 1  K0  —  1 


9i(t) 
h0(x  i) 


dt 


+ 

+ 

+ 

+ 

+ 

+ 


('d  2 n  (t  —  Xi  cos#)[(f  —  X\  cos#)2  —  (— X\  sin#)2]  .  2 
,  1  +  k  [(£  —  X\  cos#)2  +  (— X\  sinf))2]2 

rd  2 n  (t  —  X\  cos  9)[(t  —  x\  cos  9)2  +  3(— X\  sin  9)2}  2 

,  1  +  k  [(t  —  sq  cos#)2  +  (— oqsinf?)2]2 

rd  2g  (— x\  sin#)[—  (t  —  X\  cos#)2  +  (— X\  sin#)2] 


sin  6  g2(t)dt 
cos2  9  g2(t)dt 


K 


[{t  —  x\  cos  6 )2  +  (—xi  sin  9) 
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rd  2/i  (— x\  sin#)[(— x\  sin#)2  +  3(f  —  x\  cos#)2]  ._2 
,  1  +  k  [(t  —  x\  cos#)2  +  (— xi  sin#)2]2 
rd  2  id  (— aqsin#)[(— aqsin#)2  —  (t  —  aqcos#)2] 


sin  9  cos  9  g2(t)dt 
sin2  9  h2{t)dt 


1  +  k  [(£  —  x\  cos#)2  +  (— x\  sin#)2] 


212 


cos2  9  h2(t)dt 


rd  2 u  (t  —  x\  cos#)[— (— x\  sin#)2  +  (t  —  x\  cos  9)2} 

,  irr — [(t-i1cM»)»  +  (UL)T  sm()C08Wt 


— a 


oo 

yy  ' 


(101) 
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5.  Solution  to  the  Crack-Inclusion  Interaction  Problem 


Referring  to  figure  10,  the  solution  to  the  interaction  between  an  elliptical  inclusion  and  an 
arbitrarily  oriented  straight  crack  is  given  by  solving  these  four  equations  102-105  in  four 
unknowns  so  that  we  have  on  the  inclusion  (a  <  X\  <  b). 


2/i  1 

(1  +  k)  n 


(102) 
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2/i  1 

(1  +  ft)  7T 


6 


(103) 


hi(t) 
t  —  x  i 


dt  —  /i0 


fei(t) 

h0{x\) 


dt 


+  /  Gxy(xi,t)g2(t)dt+  I  Hxv(x1,t)h2(t)dt  =  -a™, 


xy\ 


xy 


and  on  the  crack  (c  <  a:2  <  d). 


Hyy(X2,  t)hi(t)dt 


2fi 


9*{t) 


(104) 


(1  +  k)  n  Jc  t-x2 


dt  =  —a™  sin2  6  —  a™  cos2  9  +  2cr™  sin  9  cos  9  , 


xy 


fb  pb 

Gl1Jx2,t)g1(t)dt+  /  H™(x2,t)hi(t)dt 


xy\ 

2/1 


rd 


_ h2(t) 

(1  +  ft)  7T  Jc  t—x2 


(105) 


dt  —  —  {a™  —  a™)  sin  9  cos  9  —  a™  (cos2  9  —  sin2  9)  , 


where  /i  and  v  are  elastic  constants  of  the  matrix;  /j0  and  vQ  are  elastic  constants  of  the  inclusion; 
k  and  k0  are  given  by  equation  20  as  a  function  of  v  or  v0,  respectively;  hQ{x i)  is  a  function  that 
describes  the  height  of  the  inclusion  as  a  function  of  X\  along  the  length  of  the  inclusion;  the 
unknown  functions  g±,  hi,  g-> ,  and  h->  are  given  by  equations  39-42;  the  Green’s  function 
transforms  G™  and  Hf-  are  given  by  equations  74-81;  and  the  Green’s  functions,  Grj  and  llrj  in 
equations  74-81  are  given  by  equations  65-69. 

5.1  Reduction  to  Well  Known  Solutions 


If  we  let  g-2  =  h2  =  0,  this  is  the  same  as  removing  the  crack  from  the  medium  so  that  we  now 
have  an  infinite  elastic  medium  containing  an  elliptical  inclusion.  With  g2  =  h2  =  0,  there  is  no 
crack  to  influence  the  inclusion  and  the  coupling  stress  terms  of  equations  102  and  103  disappear, 
and  equations  104  and  105  become  zero.  Then  we  can  remove  the  x2,  y2  coordinate  system  from 
figure  10  and  let  the  Xi,yi  coordinate  system  be  simply  the  x,  y  coordinate  system  and,  as  in 
figure  11,  place  the  origin  at  the  center  of  the  inclusion  and  normalize  so  that  on  the  inclusion 
(— 1  <  x  <  1),  we  have  two  equations  in  two  unknowns 


2/X  1  f1  gi(t)  Kq  +  1 

(1  +  k)  7T  J_x  t  —  x  ^°n0-l 


9i(t) 

h0(x) 


dt  = 


-cr. 


yy 


(106) 
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Figure  11.  Example  of  crack-inclusion  interaction:  composites, 
reinforcing  fiber  in  an  elastic  matrix.  Adapted  from 
Savalia  (22). 


2/i  1 

(1  +  K)  7T 


hi{t) 
t  —  x 


dt  —  /iQ 


hi(t) 

hQ(x) 


dt  =  —a 


OO 

xy  * 


We  can  examine  this  solution  further  for  an  elliptical  cross  section  inclusion.  We  define  the 
inclusion  height  function  hQ{x )  as 


h0{x )  =  c0V  1  -  x 2  , 

where  c0  is  some  scalar  coefficient  that  describes  the  magnitude  of  the  inclusion. 
Then  the  solution  becomes 


2/i  1  ['  9iW  ,,  _  ^o  +  l 

(1  +  k)  7T  f  -  X  ^°K0  —  1 


gig) 

c0\/l  —  X2 


dt  =  — cr 


OO 

yy  ’ 


2/r  1  f1  hi(t) 


(1  +  k)  7T  J_i  t  —  x 


,  f.  r  mo  Too 


If  we  simplify  these  expressions  by  filtering  out  the  constants  such  that 


(107) 


(108) 


(109) 

(HO) 
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(Ill) 


Cl 


(«  +  !) 

2/i 


c2 


/io(^  +  l)(«o  +  1) 
2 iicq{k0  -  1) 


(112) 


/z0(fi;  +  1) 

2/ic0 


(113) 


then  we  can  write 


1 

7T 


01  (g 

t  —  X 


dt  —  c2 


gig) 
\/l  —  a:2 


dt  =  - cicr “  . 


1 

7T 


Mg 

t  —  X 


dt  —  C3 


Mg 

\/l  —  x2 


dt  =  -ciM  • 


Then  the  solutions  to  equations  1 14  and  115  are  found  to  be 


0i(g 


Ci  O', 


yy 


1  +  c2 


t 

MM2 


Mg 


giM_ 

1  +  C3 


t 

MM 


(-1  <  f  <  1) . 


(114) 


(115) 


(116) 


(117) 


For  known  elastic  constants  and  known  loading  at  infinity,  equations  116  and  117  may  be  used  to 
determine  the  normal  and  shear  stress  fields  caused  by  the  presence  of  an  elliptical  inclusion  in  an 
infinite  elastic  medium. 


There  are  two  noteworthy  special  cases  discussed  by  Erdogan  and  Wei  (27).  First,  if  the  crack  is 
removed  as  just  discussed  but  the  stiffness  of  the  remaining  inclusion  is  zero  (i.e.,  fiQ  =  0),  then 
we  have  mathematically  the  same  situation  as  a  crack.  That  is,  an  inclusion  in  an  infinite  elastic 
matrix,  if  the  inclusion  has  zero  stiffness,  has  the  same  boundary  conditions  as  a  crack  in  an 
infinite  elastic  matrix.  This  situation  is  also  representative  of  a  void  in  an  infinite  elastic  matrix. 
In  the  case  of  an  inclusion  with  zero  stiffness,  the  solution  reduces  to  the  solution  for  a  crack  in 
linear  elastic  fracture  mechanics,  which  is  well  known  and  may  be  found,  for  example,  in  a 
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fracture  mechanics  textbook  (13).  Second,  if  the  crack  is  not  removed  but  the  stiffness  of  the 
inclusion  is  zero  (i.e.,  fi0  =  0),  then  as  before  the  boundary  conditions  in  the  zero  stiffness 
inclusion  are  the  same  as  in  a  crack:  that  is,  the  stresses  within  the  inclusion  become  zero.  Then 
the  terms  in  equations  102  and  103  derived  from  equations  56  and  57,  which  are  the  stresses  in 
the  inclusion,  become  zero.  That  is,  the  g0  J)' 1  terms  of  equations  102  and  103  become  zero,  and 
the  resulting  solution  given  by  equations  102-105  reduces  to  the  solution  for  two  arbitrarily 
oriented  cracks. 

5.2  Solution  by  Numerical  Methods 

Various  authors  have  described  solution  of  problems  in  fracture  mechanics  using  the  boundary 
integral  equation  method,  for  example,  an  early  treatment  by  Tan  and  Fenner  (23).  These 
methods  have  also  been  extended  to  solution  of  the  crack-inclusion  interaction  problem,  for 
example,  by  Tan  et  al.  (24).  This  report  does  not  provide  a  detailed  examination  of  the  solution  to 
such  problems  by  numerical  methods,  but  this  section,  adapted  from  Santare  (25),  is  included  as  a 
very  brief  introduction  to  the  solution  method  for  the  integral  equations  discussed  above. 

If  we  have  a  cracked  elastic  body  and  we  take  the  origin  of  the  coordinate  axes  at  the  center  of  the 
crack,  as  in  figure  12,  and  then  if  we  normalize  by  dividing  all  lengths  by  the  half-crack  length,  a, 
we  have  something  similar  to  the  following  Fredholm  integral  equation  of  the  second  kind: 

f  — —  dt  +  f  K(x,t)B(t)dt  =  f(x)  (118) 

J -i  t  —  x  J _i 

Observe  that  equation  118  is  similar  in  form  to  the  integral  equations  102-105  discussed  earlier. 
Here  B(t)  represents  unknown  functions  similar  to  gi(t),  h\(t),  g2(t),  and  h2(t).  K(x,  t)  is  the 
kernel  of  the  integral  equation,  the  Green’s  functions  for  the  problem  similar  to  GtJ  and  Hij, 
which  will  be  different  for  each  problem  depending  on  geometry  (singularity)  and  boundary 
conditions.  f(x)  is  some  function  dependent  on  boundary  conditions,  for  example,  in 
equations  102-105  the  tractions  applied  at  infinity  and  the  stresses  in  the  inclusion  (which  in 
equations  102-105  have  been  moved  to  the  left  hand  side  of  the  equation). 

As  in  figure  12,  we  have  a  crack,  which  we  know  must  contain  a  singularity.  So  we  know  the 
unknown  function  B(t)  must  contain  a  singularity.  To  make  the  unknown  function  easier  to 
determine,  we  can  take  out  the  singular  portion.  If  we  define  a  function  <p(t ),  which  is  bounded 
(i.e.,  not  singular)  and  well  behaved  for  —  1  <  t  <  1,  then  we  have 
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Figure  12.  Cracked  elastic  body. 


(119) 


We  know  equation  119  has  the  necessary  square  root  singularity  of  a  crack  tip  at  the  crack  tips, 
when  t  =  1  and  when  t  —  —  1.  Then  by  removing  this  singularity  from  the  unknown  function,  it 
is  easier  to  find  0(f)  than  it  would  have  been  to  find  B(t). 

To  find  0(f),  we  begin  by  discretizing  the  crack  into  evenly  spaced  finite  elements  (or  points  in 
the  one-dimensional  case  shown  here),  as  in  figure  13.  The  points  Uk+l  are  separated  by  a 
distance  h.  The  points  are  then  subdivided  by  points  0(f)  in  between  the  t2k+i  points.  Now  we 
can  use  a  numerical  interpolation  method  such  as  the  Lagrangian  interpolation  function  to 
interpolate  the  values  of  0(f)  between  the  t2k+i  points. 

Lagrange  interpolation  uses  parabolas  to  interpolate  the  curved  line  between  the  points  as  seen  in 
figure  13.  The  first  term  of  the  Lagrangian  interpolation  function  is 
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Figure  13.  Discretized  crack. 


i(.\  if,  \  [(0  ^2 k)  (I  ^2fc+l)] 

<p{t)  «  (p{t2k-l) - 2^2 -  (12°) 

Adding  more  terms  provides  a  more  accurate  interpolation.  The  second-order  accurate  Lagrange 
interpolation  approximation  then  is  given  by 


cj)(t)  = 


+ 


^  [(0  -  hk)  ( t  -  t2k+i)} 

0(f2fc— i)  2h2 

(121) 

if,  \  [(0  ~  hk-i)  (t  -  t2k+i)] 

nhk)  2h2 

(122) 

(123) 

Presumably  we  know  the  value  of  the  function  at  the  boundaries  (thus,  the  name,  boundary 
integral),  so  we  can  take  the  known  value  of  the  unknown  function  0  at  the  first  f  point  (e.g.,  at 
one  boundary)  and  the  value  of  0  at  the  last  f  point  and  interpolate  for  some  number  of  f  points 
between.  Then  take  the  interpolated  values  of  0(f)  at  known  values  of  f  and  insert  each  0(f)  and 
corresponding  f  into  equation  1 19  to  find  the  value  of  B(t)  at  that  point  f.  Then  for  equation  119 
we  have  B(t )  in  terms  of  a  set  of  unknown  numerical  values,  0(f«),  instead  of  a  set  of  unknown 
functions  0(f).  We  can  then  plug  these  values  of  B(t)  into  equation  118,  the  integral  equation, 
then  integrate.  Integrating  equation  1 18  we  get  an  algebraic  function  f{x)  in  terms  of  a  set  of 
unknown  numbers  instead  of  an  integral  equation  in  terms  of  unknown  functions. 
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Then  referring  to  figure  13,  the  solution  is 


27V+1 

E  +  VMj)K(^,ti))  m  =  /(£,)  ,  3  =  1,  2,  •  ■  •  27V  (124) 

i=l 

where  IT7*  and  Vt  are  weight  functions  determined  from  the  assumed  form  of  B(x). 


6.  Conclusions 


6.1  Stress  Intensity  Factor 

According  to  Erdogan  et  al.  (19),  when  the  stress  transformations  are  made  before  solving  the 
integral  equations,  as  in  this  report,  the  stress  intensity  factors  may  be  determined  as  follows. 
Taking  the  standard  definition  of  Modes  I  and  II  stress  intensity  factors  given  by  equations  125 
and  132,  and  if  we  have  solved  the  integral  equations  102-105  for  the  displacement  derivatives 
gi,  hi,  r/2,  and  h2  defined  by  equations  39-42,  then  we  can  determine  the  stress  intensity  factors 
using  equations  133-140.  Others,  for  example,  Erdogan  et  al.  (19)  and  Erdogan  and  Wei  in  (21), 
have  solved  the  integral  equations  using  the  technique  described  by  Erdogan  (26),  and  from  their 
results,  they  have  determined  the  stress  intensity  factors  for  crack-inclusion  interactions  of 
various  geometries  and  stiffness  ratios  between  the  shear  modulus  of  the  matrix  (ji)  and  the  shear 
modulus  of  the  inclusion  (/i0).  Their  stress  intensity  factor  results  are  tabulated  and  reported 
elsewhere  (19,  21),  so  only  qualitative  results  are  discussed  in  this  report. 

Thus  the  stress  intensity  factors  are  given  by  the  following  definitions.  At  the  a  and  b  tips  of  the 
inclusion  (a  <  x\  <  b)  in  figures  9  and  10,  the  stress  intensity  factors  are  given  by 


Ki(a)  =  lim  \J2 (a  -  aq)  a1  Jx i,0)  ,  (125) 

x\— >a 

Kn(a)  =  lim  \/2 (a  -  aq)  a1  (x i,  0)  ,  (126) 

x\— >a  y 

A7(6)  =  lim  y/2(b  ~  xi)  cr1  Ax i,  0)  ,  (127) 

xi— tb 
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Kn{b)  =  lim  a/2 {b  -  xi)  a1  (x i,0)  . 

xi—tb 


(128) 


At  the  c  and  d  tips  of  the  crack  (c  <  x2  <  d)  in  figures  9  and  10,  the  stress  intensity  factors  are 
given  by 


Ki{c)  =  lim  a/2 (c  -  x2)  crl  {x 2,0)  ,  (129) 

X2  — >C 

Kn(c)  =  lim  \/2(c~x2)  a2x Jx2,0)  ,  (130) 

X2  — >C  y 

Ki(d)  =  lim  y/2(d-  x2)  a2 Jx2,0)  ,  (131) 

X2^fd 

Kn(d)  —  lim  a/2(g?  -  x2)  /L(x2,  0)  .  (132) 


Then  the  stresses,  c/-  (i,j-x,  y,  k=  1,  2),  needed  to  determine  the  stress  intensity  factors  in 
equations  125-132  may  be  found  by  solving  the  integral  equations  102-105  for  the  displacement 
derivatives  gi,  hi,  g2,  and  h2  defined  by  equations  39-42.  Then  recognizing  that  the 
displacement  derivatives  (g,,  hu  g2,  and  h2)  along  with  the  elastic  constants  (/i  and  k)  give  the 
stress,  we  can  substitute  these  into  equations  125-132,  then  the  stress  intensity  factors  are  given 
by  the  following  equations  133-140. 

At  the  a  and  b  tips  of  the  inclusion  (a  <  x i  <  b )  in  figures  9  and  10,  the  stress  intensity  factors  are 
given  by 


KT{a)  =  lim  y/2(x[  -  a)  cq(aq)  , 


1  +  KH 


(133) 


Kn(a )  =  lim  s/2(xi  -  a)  /ii(xi)  , 

i  +  K  %i~+a 


(134) 


Kr(b)  =  lim  a/2 (b  -  Xi)  g/x/  , 


1  +  K  xi— tb 


(135) 
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Kn{b) 


T~j~— ^  lim  y/2(b-x1)  hi(xi)  . 

1  T  ^  >-6 


(136) 


At  the  c  and  d  tips  of  the  crack  (c  <  x2  <  d)  in  figures  9  and  10,  the  stress  intensity  factors  are 
given  by 


Ki{c)  =  2/i  lim  a/2(x2  -  c)  g2(a;2)  ,  (137) 

I  +  K  X2 ->c 

K  11(c)  =  lim  \/2(x2  -  c)  h2(x2)  ,  (138) 

1  +  ft  X2^fC 

KI(d)  =  -  2/'  lim  y/2(d  -  x2)  #2(2:2)  ,  (139) 

I  +  ft 

Kn(d)  =  —  lim  y/2(d  -  a:2)  h2(x2)  .  (140) 

1  +  ft  a;2^rf 


6.2  Summary 

Generally  for  a  matrix  with  a  crack  tip  near  an  inclusion,  the  crack  will  propagate  toward  the 
inclusion  if  the  inclusion  stiffness  is  less  than  the  matrix  stiffness  (including  zero  stiffness  for  a 
void  or  another  crack).  Conversely,  the  crack  will  tend  to  propagate  away  from  the  inclusion  if 
the  inclusion  stiffness  is  greater  than  the  matrix  stiffness.  Figure  14  illustrates  qualitatively  how  a 
crack  will  interact  with  an  inclusion  depending  on  the  relative  stiffness. 

We  found  that  the  stress  field  solution  for  an  edge  dislocation  has  the  properties  of  a  delta  function 
and  so  may  be  used  as  a  Green’s  function.  We  also  found  that  a  crack  may  be  represented  as  a 
continuous  distribution  of  dislocations  and  if  we  know  the  Green’s  function,  we  can  determine  the 
stress  field  for  the  crack.  So  we  used  the  stress  field  solution  for  an  edge  dislocation  as  the 
Green’s  function  for  a  continuous  distribution  of  dislocations,  which  models  a  crack  in  an  infinite 
elastic  medium.  We  also  found  that  we  can  determine  the  stress  field  in  and  around  an  elastic 
inclusion  in  an  infinite  elastic  medium  if  we  have  the  Green’s  function.  Then  also  we  can  use  the 
same  edge  dislocation  stress  field  solution  as  the  Green’s  function  to  determine  the  stresses  due  to 
the  presence  of  an  inclusion.  Using  the  Green’s  function  for  an  edge  dislocation  as  the  kernel  in  a 
set  of  integral  equations  that  describe  the  stress  field  in  and  around  an  inclusion  and  a  crack,  we 
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inclusion  stiffness  >  matrix  stiffness 


Figure  14.  The  interaction  between  a  crack  and  an  inclusion.  Schematic 
inspired  by  Pais  (27). 


can  see  from  the  form  of  a  Fredholm  integral  equation  that  we  have  a  set  of  unknown  functions 
that  must  be  identified  to  solve  for  the  stress  field  in  an  infinite  elastic  medium  due  to  the  presence 
of  a  crack  and  an  inclusion.  We  found  that  we  can  use  functions  that  describe  continuous 
distributions  of  dislocations  to  represent  a  crack  and  an  inclusion.  Then  we  use  these  unknown 
distribution  functions  as  our  unknown  functions  in  our  integral  equations.  Defining  the  form  of 
these  unknown  functions  as  the  derivatives  of  the  discontinuities  in  displacement  due  to  the 
presence  of  a  crack  in  an  elastic  medium  and  the  discontinuities  in  displacement  due  to  the 
presence  of  an  inclusion  in  an  elastic  medium,  we  can  assemble  the  integral  equations.  Finally,  if 
we  remove  the  singular  portion  of  the  unknown  functions,  to  make  them  well  behaved,  we  can 
discretize  the  crack  into  finite  elements  and  interpolate  between  the  known  values  of  the  unknown 
functions  at  the  boundaries  of  the  crack  and  inclusion  and  plug  the  results  into  the  integral 
equations,  which  we  then  solve  simultaneously  to  determine  numerical  values  for  the  unknown 
functions  in  the  integral  equations.  Then  we  have  numerical  values  for  the  stress  distribution  in 
and  around  a  crack  and  an  inclusion  and  we  can  use  these  values  to  determine  the  stress  intensity 
factors  that  result  in  an  elastic  medium  due  to  the  interaction  between  a  crack  and  an  inclusion. 
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